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collisions at s/s = 1.96 TeV 
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We present a model-independent measurement of t-channel electroweak production of single top 
quarks in pp collisions at ^fs = 1.96 TeV. Using 5.4 fb _1 of integrated luminosity collected by the 
DO detector at the Fermilab Tevatron Collider, and selecting events containing an isolated electron 
or muon, missing transverse energy and one or two jets originating from the fragmentation of b 
quarks, we measure a cross section o~(pp — > tqb + X) = 2.90 ± 0.59 (stat + syst) pb for a top quark 
mass of 172.5 GeV. The probability of the background to fluctuate and produce a signal as large as 
the one observed is 1.6 x 10 -8 , corresponding to a significance of 5.5 standard deviations. 

PACS numbers: 14.65.Ha; 12.15.Ji; 13.85.Qk 



Within the standard model (SM), top quarks are 
produced through two different types of interactions at 
hadron colliders. The largest yield is from ti pairs, 
produced via the strong interaction Pair production 
was used by the DO and CDF collaborations to establish 
the existence of the top quark in 1995 @, The second 
source of top quarks is the electroweak production of 
single top quarks, observed by the DO and CDF collab- 
orations in 2009 0, [B[ . The single top quark production 
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rate is about half that for tt pairs and the signal-to- 
background ratio is worse, so the observation of this 
process is more difficult. 

At the Tevatron, single top quarks are produced 
through W boson exchange, and accompanied by a b 
quark in the s channel (tb + ib = tb) [6(, or by both 
a b and a light quark in the t-channel (tqb + tqb = 
tqb) @, HI- A third process, usually called "associated 
production," in which the top quark is produced together 
with a W boson, has negligible cross section at the 
Tevatron 0). Prior DO and CDF publications H i] 
measured the production cross section for the sum of 
tb and tqb while assuming the relative rate between the 
two processes as predicted by the SM. This relative rate, 
however, departs from the SM prediction in several new 
physics scenarios, for example, in models with additional 
quark generations [13], new heavy bosons flavor- 
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changing neutral currents (FCNC) [l2j, or anomalous 
top quark couplings [l3l - [l5j . It is therefore of interest 
to measure the tb and tqb production cross sections 
separately, using a methodology that is independent of 
theoretical assumptions about their relative rates [HI ]. 
Furthermore, measuring the production rates of the two 
single top quark production modes to a precision compa- 
rable to that of the tt pair cross section not only tests 
the validity of the parton distribution functions of the 
proton and antiproton, but also the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements [16[ and higher-order 
corrections from quantum chromodynamics. These high- 
precision comparisons are important to check the self- 
consistency of the standard model (SM). 

This Letter builds on the techniques developed by 
the DO collaboration for prior measurements of single 
top quark production 0, 17-111, an d analyzes a larger 
dataset corresponding to an integrated luminosity of 
5.4 fb _1 . Using improved multivariate analysis (MVA) 
techniques, we report a new measurement of the tqb 
production rate that makes no assumptions about the 
tb production rate, and is therefore sensitive to various 
new physics scenarios. 

The data for this analysis were collected with a logical 
OR of many trigger conditions that results in a fully 
efficient trigger selection for the single top quark signal. 
Electron, muon, missing transverse energy ($r), and jet 
reconstruction algorithms used in this analysis have been 
described elsewhere [3, H3|- Events are selected that 
contain one jet with transverse momentum px > 25 GeV 
and at least a second jet with pr > 15 GeV, both within 
pseudorapidity \rj\ < 3.4. Events with 2, 3 or 4 jets are 
used in the analysis. Events are also required to contain 
exactly one isolated high-py electron or muon that origi- 
nates from the pp interaction vertex and satisfies the 
following acceptance criteria: for the electron \rj\ < 1.1 
and pr > 15(20) GeV for events with 2 (3 or 4) jets; for 
the muon |?7| < 2.0 and pr > 15 GeV. The |5y is required 
to be in the range of (20, 200) GeV for events with 2 jets 
and (25, 200) GeV for events with 3 or 4 jets. 

The main background contributions arise from 
W bosons produced in association with jets, tt pairs, and 
multijct production in which a jet with high electromag- 
netic content mimics an electron, or a muon contained 
within a jet originating from the decay of a heavy-flavor 
quark (b or c quark) appears isolated. Diboson [WW, 
WZ, ZZ) and Z+jets processes add small additional 
contributions to the background. To limit the number 
of events originating from multijet production that enter 
our candidate samples through ^?y mismeasurements, we 
require the $t to be not aligned with or opposing the 
direction of the lepton or the leading jet and require 
high total scalar transverse energy of the lepton, jets 
and TpT (Ht)- Additional requirements are placed on the 
transverse energy of all jets in all events and the second 
leading jet py in two-jet events to suppress soft-scattering 



background. 

To enhance the fraction of signal events in data, 
one or two of the jets in each event are required to 
originate from the hadronization of long-lived b hadrons, 
as inferred from a neural network 6-tagging algorithm 
("b-tagged jet") [2l|. The data are divided into 
six mutually-exclusive subsamples to take advantage 
of different signal-to-background ratios and different 
sources of background. This separation is based on jet 
multiplicity (2, 3, or 4 jets), and on the number of fa- 
tagged jets (1 or 2 6-tags). The signal acceptance of 
all selections corresponds to (2.0 ± 0.3)% for tqb and 
(2.9 ± 0.4)% for tb production, as estimated using Monte 
Carlo (MC) simulations. 

Single top quark signal events are modeled using the 
COMPHEP-based MC event generator singletop [22I ]. 
The kinematics of the generated events closely match 
those predicted by next-to-leading-order (NLO) calcula- 
tions [23[ . pythia [24[ is used to model the hadronization 
of all generated partons. The tt, W+jets, and Z+jets 
backgrounds are simulated using the ALPGEN leading-log 
MC event generator [2^], with pythia used to model 
hadronization. The tt background is normalized to 
the predicted cross section [l[ for a top quark mass of 
172.5 GeV. The normalization of the W/Z backgrounds 
is obtained by scaling the ALPGEN cross sections to 
account for NLO effects [j^ . Diboson backgrounds are 
modeled using pythia. 

The presence of additional pp interactions is modeled 
by events selected from random beam crossings matching 
the instantaneous luminosity profile in the data. All MC 
events are passed through a GEANT-based simulation (2?| 
of the DO detector. Small corrections, derived using 
samples depleted of single top quark events, are applied 
to all reconstructed objects to improve the agreement 
between data and simulation. In particular, any mismod- 
eling of the pseudorapidity of jets or of the angular 
separation between the two leading jets in the W+jets 
sample is adjusted by rcweighting the MC. 

The multijet background is modeled using data 
containing non-isolated leptons. In the electron channel, 
the transverse momentum distribution of the lepton 
is reweighted to match the distribution in background 
events that pass the candidate-selection criteria. To 
increase the statistics in the muon channel, the jet closest 
to the muon is removed and the $t recalculated. The 
overall normalization of the W^+jets and the multijet 
backgrounds is obtained by comparing the background 
expectation to data in three sensitive variables: lepton 
transverse momentum, ^?y, and the transverse mass 
of the W boson. This normalization is found after 
subtracting the contributions from the small backgrounds 
{tt, Z+jets, and dibosons), before 6-tagging, when the 
expected signal-to-background ratio (S:B) for the entire 
sample is S:B 1:279. 

The probability for the 6-tagging algorithm to identify 
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a jet as originating from a b quark is measured in a 
sample of data with jets containing muons and compared 
to the corresponding efficiency measured in a MC sample. 
From this comparison we derive correction factors that 
we parametrize as a function of jet flavor, pt, and 77 and 
apply to all the MC samples used in this analysis. After b- 
tagging, we check the normalization of the Wbb and Wcc 
samples in a sample that has no overlap with the one 
used in the single top quark cross section measurement 
and find the normalization to be consistent with unity. 
We assign a 12% uncertainty on the normalization of the 
Wbb and Wcc samples to account for variations in the 
Wcj (j — u, d, s) cross section and the Wbb to Wcc cross 
section ratio used in this study. 

We consider systematic uncertainties that arise from 
the theoretical and detector modeling of the background 
processes and from the subsequent corrections applied to 
that model. Most uncertainties affect only the normal- 
ization of the distributions, but the corrections for the jet 
energy scale (JES), the jet energy resolution (JER), and 
the reweighting of the kinematic distributions in IT±jcts 
events also modify the differential distributions. The 
largest uncertainties come from the JER, corrections to 
the 6-tagging efficiency, and the corrections for the jet- 
flavor composition in M^+jets events, with smaller contri- 
bution from JES, MC statistics, integrated luminosity, 
and trigger uncertainties. The total systematic uncer- 
tainty on the background is 11%. 

We select 8,471 events and expect 238 ± 28 tqb events. 
This corresponds to a S:B of approximately 1:33 for the 
sample with 1 6-tag and 1:50 for the sample with 2 b- 
tags. Table Q] lists the number of events expected and 
observed for each process as a function of jet multi- 
plicity. We also define two control samples to check 
the background model components separately for the 
two main backgrounds: IF±jets and ti. The W+jets- 
dominated sample has low Ht, exactly two jets, and 
only one 6-tagged jet. The tt-dominated sample has high 
Ht, exactly four jets, and one or two ^-tagged jets. We 
find good agreement for both normalization and shape in 
most of distributions studied for the cross-check samples 
and also for the candidate sample before 6-tagging, when 
the signal contribution is negligible. 

We further improve the discrimination between signal 
and background by employing MVA techniques as 
described in 18j. We use two of our previous methods: 
boosted decision trees (BDT) [28| and Bayesian neural 
networks (BNN) j2^|. In addition, we use a third method, 
neuroevolution of augmented topologies (NEAT) [30| . 
All three methods use the same data and background 
model, consider the same sources of systematic uncer- 
tainty, and rely on variables with signal-background 
discrimination power chosen from a common set of 
well modeled variables from five categories [HI I32I ]: 
single object kinematics, global event kinematics, jet 
reconstruction, top quark reconstruction, and angular 



TABLE I: Number of expected and observed events in 
5.4 fb _1 of integrated luminosity, with uncertainties including 
both statistical and systematic components. The tb and tqb 
contributions are normalized to their SM expectations. 



Source 



2 jets 



3 jets 



4 jets 



tb 
tqb 

it 

W+jets 



104 ± 16 44 ± 7.8 13 ± 3.5 

140 ± 13 72 ± 9.4 26 ± 6.4 

433 ± 87 830 ± 133 860 ± 163 

3,560 ± 354 1,099 ± 169 284 ± 76 



Z±jets & dibosons 400 ± 55 142 ±41 35 ± 18 
Multijets 277 ± 34 130 ± 17 43 ± 5.2 

Total prediction 4,914 ± 558 2,317 ± 377 1,261 ± 272 
Data 4,881 2,307 1,283 



correlations. The BDT method uses a common set 
of 50 variables for all analysis channels, while the 
NEAT method uses the TMVA implementation of the 
"RuleFit" algorithm [H, SH to select the 30 most 
sensitive kinematic variables in each channel. New for 
this analysis, the BNN method constructs a discriminant 
using the lepton and jet 4-vectors, the J^t 2-vcctor, and 
variables that include lepton charge and 6-tagging infor- 
mation, for a total of 14, 18 and 22 variables for events 
with 2, 3 and 4 jets. 
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FIG. 1: Distribution for tb and tqb processes for (a) trans- 
verse momentum of the leading 6-tagged jet, (b) centrality, 
defined as the ratio of the scalar sum of the pr of the jets 
to the scalar sum of the energy of the jets in the event, (c) 
lepton charge multiplied by the leading non-6-tagged quark 
jet pseudorapidity and (d) total transverse energy. All the 
distributions are normalized to unity. 

Each method is optimized to maximize the sensi- 
tivity to the tqb signal by treating the tb process as a 
background component with normalization given by the 
SM cross section Q . Figure Q] shows the distributions of 
some of the variables that provide the greatest discrim- 
ination between the tqb and tb processes. To further 
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improve the sensitivity, the analyses are performed 
separately on the six mutually exclusive subsamplcs 
defined before. For all three methods the output 
variables are saved in histograms with a binning chosen 
to ensure that there are enough events to limit the uncer- 
tainties due to MC statistics. 

Even though the three MVA techniques use the same 
data sample, they are only 70% correlated with 
each other. We therefore combine these methods using 
an additional BNN algorithm (BNNComb) that takes 
as input the individual output discriminants of the 
BDT, BNN, and NEAT methods, and produces a single 
combined output discriminant. Figure [2] shows compar- 
isons between the t-channel signal, the background 
model, and data for the combined discriminant, which 
leads to a more precise measurement of the cross section. 
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FIG. 2: Comparison of the signal and background models 
to data for the combined t-channel discriminant for (a) the 
entire discriminant range and (b) the signal region. The bins 
have been ordered by their expected S:B. The single top quark 
contributions are normalized to the measured cross sections. 
The t-channel contribution is visible above the hatched bands 
that show the uncertainty on the background prediction. 

The single top quark production cross section is 
measured using a Bayesian approach as in 0, 17, [lH ]. 
We follow the approach of [19( and construct a two- 
dimensional (2D) posterior probability density as a 
function of the cross sections for the tqb and tb processes. 
The output discriminants for the signals, backgrounds, 
and data are used to form a binned likelihood as a 
product over all six analysis channels and all bins. No 
constraint is imposed on the relative rates of tb and 
tqb production. We assume a Poisson distribution for 
the observed number of events and uniform prior proba- 
bilities with positive values for the two signal cross 
sections. We integrate over the systematic uncertainties 
which are described by Gaussian priors that preserve 
the correlations between bins and channels. The tqb 
cross section is then extracted from a one-dimensional 
posterior probability density obtained from this 2D 
posterior by integrating over the tb axis, thus not making 
any assumptions about the value of the s-channel cross 
section. Similarly, the tb cross section is obtained by 
integrating over the tqb axis. Ensembles of datasets 
generated at several different cross section values are 
used to verify the linearity of the cross section extraction 



procedure. 

Figure [3] shows the 2D posterior probability density 
for the combined discriminant together with predictions 
from the SM || and various beyond-the-SM scenarios: 
four-quark-gencrations with CKM matrix element |Vj s | = 
0.2 [lOf, top- flavor model with new heavy bosons at a 
scale m x = 1 TeV and FCNC with an up-quark/top- 
quark/gluon coupling k u /A = 0.036 fl2l ]. 
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FIG. 3: Posterior probability density for tqb vs tb single top 
quark production in contours of equal probability density. 
The measured cross section and various theoretical predic- 
tions are also shown. 

The measured cross sections of a(pp —> tqb + X ) = 
2.90 ± 0.59 pb and a{pp -+tb + X) = 0.98 ± 0.63 pb are 
in good agreement with the SM expectation for a top 
quark mass of 172.5 GeV @. The uncertainty includes 
both statistical and systematic sources. The cross section 
for t-channel single top quark production is the most 
precise measurement of an individual single top quark 
production channel to date with an uncertainty of 20%. 

The significance of the t-channel cross section 
measurement is computed using a log-likelihood ratio 
approach [H, [p| which tests the compatibility of the 
data with two hypotheses: a null hypothesis where 
there is only background and a background plus signal 
hypothesis, where the number of signal events corre- 
sponds to the theoretical cross section. New for this 
analysis is the computation of the distributions for these 
two hypotheses given by an asymptotic Gaussian approx- 
imation [3o| . With this approximation we compute 
for the first time, the significance of the measured tqb 
cross section independently of any assumption on the 
production rate of tb. We estimate the probability of the 
background to fluctuate and produce a signal as large 
as the one observed to be 1.6 x 10 -8 , corresponding 
to a significance of 5.5 standard deviations (SD). The 
expected significance is 4.6 SD. 

The presence of the t-channel signal is visible in 
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Fig. |U which shows comparisons between the data, the 
background model, and the t-channel signal for a sample 
selected with S:B > 0.32 using the combined discrim- 
inant. Four kinematic variables with large discrimi- 
nating power for t-channel single top quark production 
are shown: Ht, the reconstructed top quark mass, the 
lcpton charge multiplied by the pseudorapidity of the 
leading non-6-tagged jet, and the t-channel spin corre- 
lation in the optimal basis [32], i.e., the cosine of the 
angle between the light quark jet and the lepton, both in 
the reconstructed top quark rest frame. 
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FIG. 4: (a) Total transverse energy, (b) reconstructed top 
quark mass, (c) lepton charge multiplied by the leading non-6- 
tagged quark jet pseudorapidity, and (d) t-channel top quark 
spin correlation (see text) for a sample selected with S:B > 
0.32 using BNNComb discriminant. The t-channel and s- 
channel contributions have been normalized to their measured 
cross sections. 

The measured cross section depends on the assumed 
mass of the top quark, m t . The dependence is studied 
by repeating the analysis on MC samples generated at 
different values of m t . Table HTl summarizes the measured 
cross section for different top quark masses. 

TABLE II: Variation of the measured cross section for tqb and 
tb single top quark production with the top quark mass. 



m t 


170 GeV 


172.5 GeV 


175 GeV 


tqb 
tb 


2.80ig;SI 
1.3H8& 


2.90t ;g 


9 rn + 0.58 

o.65±g : gs 



In summary, we have presented a model-independent 
measurement of the t-channel single top quark 
production cross section using 5.4 fb -1 of integrated 
luminosity recorded by the DO experiment at the 
Fcrmilab Tevatron accelerator. For a top quark mass 
of m t = 172.5 GeV, the measured cross section is 
a(pp —¥ tqb + X) = 2.90 ± 0.59 pb. This is the most 
precise measurement to date, and is in good agreement 



with SM expectations. 
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